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Abstract 

The 3 He( 3 He, 2p) 4 He and 3 H( 3 H, 2n) 4 He reactions are studied in a micro- 
scopic cluster model. We search for resonances in the 3 He+ 3 He and 4 He+p+p 
channels using methods that treat the two- and three-body resonance asymp- 
totics correctly. Our results show that the existence of a low-energy resonance 
or virtual state, which could influence the 7 Be and 8 B solar neutrino fluxes, 
is rather unlikely. Our calculated 3 He( 3 He, 2p) 4 He and 3 H( 3 H, 2ra) 4 He cross 
sections are in a good general agreement with the experimental data. 
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I. INTRODUCTION 



One way to test models of the solar interior is by observing the neutrinos generated by 
the nuclear reaction network which is the solar energy source [1]. As a striking and exciting 
fact, all terrestrial solar neutrino experiments observe fewer neutrinos than predicted by 
standard solar models [2]. The picture which emerges from the various experiments is [3] 
that the 8 B flux is about half of its predicted value (0§ = O.40f s ' M ), while the 7 Be neutrinos 
appear to be completely missing (0 7 = 0). Here SSM refers to the standard solar model 
of Bahcall [2]. The 7 Be nuclei, which are the seeds of both the 7 Be and 8 B neutrinos, are 
produced in the 4 He( 3 He, 7) 7 Be reaction. This reaction competes with the 3 He( 3 He, 2p) 4 He 
process, which is the final step of the first branch of solar hydrogen burning (ppl chain) 
[4]. If the cross section of the 3 He( 3 He, 2p) 4 He reaction were larger than believed then the 
probability of the 4 He( 3 He, 7) 7 Be branch would be smaller, and hence the 08 and 07 fluxes 
would be suppressed (without, however, significantly affecting the unexpected 07/08 ratio 
deduced from the neutrino experiments [5]). To increase the 3 He( 3 He, 2p) 4 He reaction rate, 
a hypothetical resonance at low energies had been suggested [6] and looked for in various 
experiments (see [7] for the most recent work). 

The importance of the 3 He( 3 He, 2p) 4 He reaction has led to continued experimental efforts 
to measure the cross section down to solar energies. In the latest experiment [7] the LUNA 
collaboration measured the cross section down to E cm = 20.76 keV, which is well within 
the region of the most effective solar energies (Gamow window). Although, they do not see 
any evidence of a possible resonance, the existence of such a state at still lower energies 
cannot be a priori ruled out yet. In the present work we study the 3 He( 3 He, 2p) 4 He and 
the mirror 3 H( 3 H, 2n) 4 He reactions in a microscopic cluster model. We search for signs of 
possible resonances and study the energy dependence of the reaction cross sections. 



II. MODEL 



The 3 He( 3 He, 2p) 4 He and 3 H( 3 H, 2n) 4 He reactions have already been studied previously 
within microscopic cluster models [8,9]. We use the same model, but with extended and 
hence more realistic model spaces. Additionally we put special emphasis on the description 
of the few-body dynamics. Our model starts out with a resonating group model (RGM) 
wave function for the { 4 He + p + p, 3 He + 3 He} coupled system. 
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Here A is the intercluster antisymmetrizer, the <3> cluster internal states are translationally 
invariant harmonic oscillator shell model states (a = 4 He and h = 3 He), the p vectors are 
the different intercluster Jacobi coordinates, l\ and l 2 are the angular momenta of the two 
relative motions, L and S are the total orbital angular momentum and spin, respectively, J is 
the total angular momentum, ix = (— = (— l)' 3 is the parity, and [...] denotes angular 
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momentum coupling. The sum over S,l±,l2, and L includes all physically relevant angular 
momentum configurations. Using (1) as a trial function in the six-nucleon Schrodinger 
equation, we arrive at an equation for the intercluster relative motion functions x- F° r the 
mirror { 4 He + n + n, 3 H + 3 H} system we use a wave function similar to (1). The harmonic 
oscillator size parameters of the internal cluster states are chosen to stabilize the cluster 
energies. We use the Minnesota (MN) effective interaction [10] in all calculations. It puts 
the threshold of the 3 He + 3 He and 3 H + 3 H thresholds at 17.11 MeV and 15.57 MeV energies, 
respectively relative to the 4 He + N + N threshold. Like in the previous studies [8,9], these 
values are significantly larger than the experimental thresholds (12.86 MeV and 11.33 MeV, 
respectively). One might expect that the incorrect reproduction of the threshold energies 
will effect the relative weight of the 4 He + N + N channels in the full wave function. We 
will briefly discuss this problem below. 

The relative motions are the most important degrees of freedom in the problem, so special 
care has to be taken in order to ensure that their dynamics are described properly. What 
makes the description of the 3 He( 3 He, 2p) 4 He and its mirror reactions difficult is the fact 
that there are three particles in the final state continuum. Currently we cannot treat the full 
three-body continuum problem properly, thus our model is only an approximate description 
of the reactions. However, the existence of a low-energy resonance in the 3 He( 3 He, 2p) 4 He 
reaction is a question which can be studied in a rigorous way. If such a resonance existed, 
it has to originate from either the 3 He + 3 He or the 4 He + p + p channels. As we will show 
below, one can investigate the existence of such a resonance in both channels separately 
and properly. Although we can treat the scattering continuum only approximately, we will 
nevertheless calculate the low-energy reaction cross sections of the 3 He( 3 He, 2p) 4 He and its 
mirror reaction. 

III. SEARCHING FOR RESONANCES IN 4 He+p+p AND 3 He+ 3 He 

In order to avoid any ambiguity in the recognition of a resonance in real-energy observ- 
ables, we search for resonances on the complex energy plane. Resonances are defined as the 
complex-energy solutions of the Schrodinger equation, which correspond to the singularities 
of the scattering matrix. Thus, we search for poles of both the 3 He + 3 He and 4 He + p + p 
S matrices. 

Interestingly, the 4 He + p + p case is easier to deal with, despite its three-body nature. 
The reason is that a low-energy resonance in the 3 He + 3 He channel corresponds to a high- 
lying narrow state in the 4 He + p + p channel, which, if it exists, can be easily identified. In 
order to be able to describe three-body resonances, we apply the complex scaling method 
(CSM). This method has already been used previously to search for resonances in 6 He, 6 Li, 
and 6 Be in an a + N + N model [11]. However, the search had been restricted to low-lying 
resonances below the 3 He+ 3 He threshold and therefore does not shed light on the problem 
at hand here. In the present search we concentrate on the high-energy region around the 
3 He + 3 He channel threshold. In the CSM the complex scaling transformation is performed 
on the original Hamiltonian. This transformation acts in coordinate space on a function 
/(r) as 

0(9) f(r) = e 3 ^ 2 /(re i0 ), (2) 
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where 9 is the parameter of the transformation. For real 9 values the U (9) transformation 
results in a rotation into the complex coordinate plane. The spectrum of the complex-scaled 
Hamiltonian 

H e = U(9)HU-\9) (3) 

is connected to the spectrum of the original H by the following theorem [12]: (i) the bound 
eigenstates of H are eigenstates of H e , for any value of 9 within < 9 < ir/2; (ii) the 
continuous spectrum of H is rotated by an angle 29; (hi) the complex generalized eigenvalues 
of He, e res = E r — iT/2 (with E r , T > 0), belong to its proper spectrum, with square- 
integrable eigenfunctions, provided 29 > |arge rcs |. These complex eigenvalues coincide with 
the S'-matrix pole positions. 

We would like to emphasize that this method treats three-body resonances in a rig- 
orous way. The only approximation we make in the present work is that we solve the 
complex-scaled Schrodinger equation on a finite basis. Namely, we assume that x i n (1) can 
be expanded in terms of products of Gaussian functions, like p l i exp[— (pi/^fi) 2 ]Yi imi (pi) ■ 
p l 2 exp[— (p2 /lj) 2 ]Yi 2 m 2 (P2) , where li and I2 are the angular momenta in the two relative mo- 
tions, respectively, and the widths 7 of the Gaussians are the parameters of the expansion. 
The expansion coefficients are determined from the projection equation (5^e\H e — e\^e) = 0. 
As the complex-scaled wave function of a resonance is square-integrable, our finite basis ap- 
proximation works for resonances as well as for bound states. 

As an illustrative example of our results, Fig. 1 shows the spectrum of the complex-scaled 
Hamiltonian for the J w = + state of 4 He + p + p. We choose a small rotation angle 9 in 
order to avoid any numerical instability. One can see that the well-known ground state 
resonance of 6 Be is reproduced. A low-energy resonance in the 3 He( 3 He, 2p) 4 He system with 
small width (which would be relevant for the astrophysical problem), if existed, would be 
situated above Ke(E) ps 10 MeV, close to the real energy axis. There is no such narrow state 
present in our model. We do not find such a state in other J n partial waves either. Our 
results agree with other calculations performed in a three-body 4 He + N + N model assuming 
structureless 4 He [13]. As Ref. [13,14] shows, all high-lying resonances in the 4 He + N + N 
systems have large widths. 

In the 3 He + 3 He channel we cannot use the complex scaling method. A very low-energy 
resonance would always be mixed up with the rotated continuum, making its unambiguous 
identification hopeless. Here we use a direct analytic continuation of the S matrix to complex 
energies [15]. We can use this method because, unlike in the case of the three-body 4 He + 
p + p system, the two-body scattering wave functions can easily be generated with the 
correct asymptotics. To calculate these scattering wave functions we use the Kohn-Hulthen 
variational method of Ref. [16]. Once again, our method treats resonances in a rigorous way. 

We have searched for low-energy narrow resonances in 3 He + 3 He and found none. We 
must note, however, that our 3 He + 3 He model may be too simple. In the J w = + state, 
which is the most likely candidate to have a narrow resonance, the singlet S'-wave N — N 
interaction has a dominant effect because the Pauli principle forces the unpaired neutrons 
inside the 3 He clusters to form a 1 S'o state between them. In the 3 He + 2 H system the role 
of the triplet forces is known to be very strong, causing the existence of the low-lying 3/2+ 
resonance [17]. The triplet forces, which have negligible effects in our present model, could 
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play a role if the small D-state admixture of 3 He were to be considered. In such a case one 
can have a contribution from coupling between the | (1,1/2)1/2,0 1/2; (1,1/2)3/2,2 1/2 j 

configurations in the 3 He+ 3 He system. Here the (Sd, S P )S, I I coupling scheme is used, with 
Sd and S p being the deuteron and proton spins, respectively, S is the total intrinsic spin, / is 
the angular momentum between the deuteron and the proton, I is the total spin of 3 He, and 
the brackets denote angular momentum coupling. Such a model would require a four-cluster 
description of 3 He + 3 He = (d + p) + (d + p), which is beyond the scope of the present work. 
The effects of the internal D states in 3 He on the 3 He + 3 He system will be studied in the 
future [18]. 

In many respects 3 He + 3 He is similar to the n+n system. We know that there is a virtual 
(antibound) state present in the n + n system, with the wave number k = —27 (7 > 0) 
and energy E = —Ey (Ey > 0) of the ^-matrix pole, respectively. Such a state results 
in a cross section which is divergent at the unphysical negative pole energy and behaves as 
a ~ 1/(E + Ey) at positive energies. The effect of such a hypothetical state in 3 He + 3 He 
on the 3 He( 3 He, 2p) 4 He cross section was discussed in Ref. [19]. It would lead to a cross 
section that increases with decreasing energy, mimicking the effect of electron screening. 
A closer look at the problem shows that pure virtual states (with pure imaginary wave 
number) cannot be present in Coulombic systems [20]. The Coulomb interaction creates 
two poles from the one virtual-state pole and moves them away from the imaginary /c-axis 
to k — ±/t — 27 (7 > 0, k < 7). Such states can still have observable effects, like in the p + p 
system, because these conjugate poles are roughly at the same distance from the physical 
energies. We have searched for such states in 3 He + 3 He, and found no unambiguous evidence 
for their existence close to the imaginary k axis. 

In summary, we do not find any S'-matrix poles in either 4 He + p + p or 3 He + 3 He that 
could cause strong observable effects in the 3 He( 3 He, 2p) 4 He cross section. For 4 He + p + p 
our result is probably the best, one can currently achieve. Further studies of 3 He + 3 He in 
a (d + p) + (d + p) model or in a full six-body dynamical model would be interesting. 



IV. THE 3 He( 3 He,2p) 4 He AND 3 H( 3 H,2n) 4 He REACTION CROSS SECTIONS 

3 He( 3 He, 2p) 4 He is the only solar nuclear reaction whose cross section has been measured 
within the solar Gamow window (around ^20 keV). At such low energies the cross section, 
measured in the laboratory, is enhanced due to screening effects by the electrons present 
in the target atoms [21]. This electron screening effect has definitely been identified in the 
latest LUNA data for 3 He( 3 He, 2p) 4 He. For applications in the solar models the electron 
screening enhancement has to be subtracted from the data. 

Currently the 3 He(<i,p) 4 He reaction is the one for which the enhancement of the low- 
energy fusion cross section due to electron screening is studied best. In agreement between 
experiment [22] and theory [23-25] it appears that the screening enhancement for this re- 
action, in which deuterons collide with an atomic 3 He gas target, is well described in the 
adiabatic limit. In this case the electron screening can be represented by a constant shift of 
the collision energy by the screening energy U e which is given by the difference of the elec- 
tronic binding energy of the united atom and the sum of the asymptotic fragments. Applied 
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to the 3 He + 3 He reaction, the screening energy in the adiabatic limit is U e = 240 eV. We 
will use this value in the following, but we note that fits to the LUNA data might indicate 
a somewhat larger screening potential [7,33]. These fits had to make assumptions about the 
energy dependence of the bare-nuclear 3 He( 3 He, 2p) 4 He S factor. This has motivated us to 
perform calculations for the bare reaction cross sections of 3 He( 3 He, 2p) 4 He and the mirror 
reaction 3 H( 3 H, 2n) 4 He. 

As we mentioned, a microscopic description of 3 He( 3 He, 2p) 4 He, which handles the full 
three-body final state rigorously, is currently not feasible. Here we use a simplified version 
of the continuum-discretized coupled channel method [26] to describe 4 He + JV + JV. In this 
method the total energy available for 4 He + N + N is divided between the 4 He — JV and 
( 4 He, JV) — JV relative motions in the ( 4 He, JV)JV configuration, and between the JV — JV and 
(JV, JV) — 4 He motions in the (JV, JV) 4 He configuration. Within the two-cluster subsystems, 
( 4 He, JV) or (JV, JV), the continuum energy is discretized, and the remaining energy appears 
as the scattering energy in the ( 4 He, JV) — JV and (JV, JV) — 4 He two-body systems. The 
system of coupled channels is built up from the ( 4 He, JV) — JV and (JV, JV) — 4 He channels 
containing the various discretized-energy states in the two-cluster subsystems. 

Generally, the discretization of the continuum in the two-body subsystems is done in 
equidistant k bins, and proper continuum states are used. Here we adopt a simpler approach. 
We discretize the continuum on a finite square-integrable basis. Thus, our discretized states 
are pseudo-bound states with square-integrable wave functions and positive energy. The 
discretization is performed by choosing Nd basis states with ranges that increase following a 
geometric progression. By varying the total range of the basis and Nd one can achieve very 
different discretization patterns, e.g., dense or sparse discretization, discretizations including 
only low energies or allowing high off-the-energy-shell states also, etc. This way one can test 
the sensitivity of the calculated reaction cross section on the continuum discretization, and 
see if this approximation is reasonable or not. We typically use Nd = 5 — 10 and solve the 
coupled-channel scattering problem by using the method of Ref. [16]. 

Figure 2 shows our results for 3 He( 3 He, 2p) 4 He and 3 H( 3 H, 2n) 4 He together with the 
available experimental data. In order to get rid of the the trivial exponential dropping of 
the cross sections caused by the Coulomb penetration, we use the astrophysical S-factor 
parametrization 



Our curves in Fig. 2 come from a continuum discretization that proved to be the most 
stable at the ( 4 He, JV) — iV and (JV, N) — 4 He two-body scattering level. We also tried other 
discretization patterns and found that the absolute normalization of the cross section curves 
depend somewhat (10-20%) on the chosen discretization, but the shapes of the curves remain 
very similar. Our results are close to those of Ref. [9], where a similar model was used. Our 
full model space is roughly 5-10 times bigger than in Ref. [9], which allows us to use much 
more flexible continuum discretizations. Nevertheless, all our results seem to be similar to 
Ref. [9], e.g., we also find that the channels involving the 4 He + JV states with J n = l/2~ 
have a dominant role. 

Fig. 3 shows the effect of electron screening on our calculated 3 He( 3 He, 2p) 4 He S factor 
with U e = 240 eV screening potential. A nice agreement is observed with the low-energy 
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LUNA data. 

The overall agreement between our results and the experimental data is considered to 
be good. We see, however, a marked disagreement with the 3 H( 3 H, 2n) 4 He data at very 
low energies. The energy dependence of our calculated S factor is different from the precise 
Los Alamos data. One possible explanation of this discrepancy could be our approximate 
treatment of the three-body final state. However, this is not supported by our finding that 
the shape of the S(E) curve is rather insensitive to the way the discretization is done. 
Nevertheless, an improved model with the full three-body treatment of the final state would 
be desirable. As we mentioned, the thresholds of the 4 He+n+n and 3 H+ 3 H channels are too 
far from each other in our model. In order to see if this may affect the energy- dependence 
of the S factor, we made some test calculations. We artificially modified the energies of 
the 3 H clusters to reproduce the experimental threshold energy difference. This changed 
the absolute normalization of the S factor (as the 3/2~ and especially the l/2~ discretized 
states moved closer to the 3 H + 3 H threshold) but not its shape. 

V. CONCLUSION 

We have studied the 3 He( 3 He, 2p) 4 He and 3 H( 3 H, 2n) 4 He reactions within the micro- 
scopic cluster model using significantly larger model spaces than previously employed. Our 
motivation and results have been twofold: 

We searched for signs of possible low-energy resonances in 3 He( 3 He, 2p) 4 He, which could 
have important effects on the 7 Be and 8 B solar neutrino fluxes. The 3 He + 3 He and 4 He + 
p + p channels were studied separately, which allowed us to use methods that can treat the 
two- and three-body asymptotics in a rigorous way. We extended the two- and three-body 
scattering matrices to complex energies and searched for their poles. We have not found 
any indication for the existence of a low-energy resonance or virtual state that could cause 
observable effects in the cross section. Thus, it is unlikely that a yet unobserved resonance 
around the threshold energy in the 3 He( 3 He,2p) 4 He reaction might affect this important 
solar reaction cross section. 

We calculated the cross sections of the 3 He( 3 He, 2p) 4 He and 3 H( 3 H, 2n) 4 He reactions 
in the continuum-discretized coupled channel approximation. Our results are in a good 
general agreement with available data, except for the very low-energy 3 H( 3 H, 2n) 4 He cross 
section, where we observe a systematic deviation from the most precise measurement. Our 
test calculations show that the energy dependence of the cross sections is hardly influenced 
by the details of the continuum discretization, but might be caused by the approximate 
treatment of the 3-body continuum. Here improvements are certainly warranted. 
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FIGURES 
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FIG. 1. Energy-eigenvalues of the complex scaled Hamiltonian of the J n = + states 
4 He + p + p. The dots are the points of the rotated discretized continuum, while the circle 
the ground state resonance of 6 Be. 
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FIG. 2. Astrophysical S factors for the (a) 3 He( 3 He, 2p) 4 He and (b) 3 H( 3 H, 2n) 4 He reactions. 
The experimental data are taken from (a) [7] (filled circle, and filled triangle), [27] (open circle), 
[28] (open triangle) and (b) [29] (filled circle), [30] (cross), [31] (square), and [32] (triangle). The 
solid curves are our results. 
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FIG. 3. The same as Fig. 2(a), except that the effect of electron screening on our low-energy 
theoretical curve is shown by the dashed line. The adiabatic screening potential, U e = 240 eV, is 
used. 
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